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^ • The structure of a quantized vortex in a Bose-Einstein Condensate is investi- 

gated using the projection method developed by Peierls, Yoccoz, and Thouless. 
This method was invented to describe the collective motion of a many-body 
system beyond the mean-field approximation. The quantum fluctuation has 
been properly built into the variational wave function, and a vortex is de- 
\ scribed by a linear combination of Feynman wave functions weighted by a 

Gaussian distribution in their center positions. In contrast to the solution 
of the Gross-Pitaevskii equation, the particle density is finite at the vortex 
axis and the vorticity is distributed in the core region. 

PACS numbers: 03.75.Fi, 67.40. Vs. 
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1. INTRODUCTION 

X . 

sation in trapped atomic gases Jil'i'ifl the generation of quantized vortices 
has become the subject of numerous theoretical and experimental studies. 
Two different schemes have been successfully implemented in magnetically 
trapped rubidium gases. The first one, which is proposed by Williams and 
Holland ,i uses a laser beam to imprint a phase winding between the two 
components in a binary mixture of condensates^ One of the components ro- 
tates around the other one with a quantized circulation. The second scheme 
creates vortices by stirring the condensate with a laser beam which alters 
the trapping potential. i In contrast to the vortices in the superfluid 4 He, 
the size of the vortex core in a gaseous condensate is orders of magnitudes 
larger, and can be observed optically. It is the structure of a singly quantized 
vortex generated by the second scheme that is studied in this work. 

In the framework of the mean-field theory developed by Bogoliubov,H 
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the ground state of a weakly interacting Bose gas can be well described by 
the Gross-Pitaevskii (GP) equationflH3 The iV-body wave function of the 
ground state is, to a good approximation, a product of N identical single- 
particle wave functions. This single-particle wave function, called the con- 
densate wave function, is the solution of the GP equation. For currently 
available experiments, the quantum depletion, which characterizes the frac- 
tion of atoms not in the condensate due to interactions, estimated using the 
local density approximation is typically about 1% or lessH 

The GP equation is often being generalized to describe the wave function 
of a collective excitation. In the case of a singly quantized vortex along the 
z axis in the center of a trapped gas, the mean-field solution corresponds to 
putting all the atoms in the single-particle state with angular momentum 
£ z = h. The condensate wave function has to vanish on the axis of the 
vortex because the phase of the condensate wave function is not well-defined 
on the axis. The modulus of the condensate wave function quickly rises from 
zero at the axis to its maximum value over a distance called healing length, 
which defines the size of the core of a vortex. However, this is the same 
length scale in which quantum fluctuations are important It is not clear 
to what extent the mean-field approximation remains valid in the vortex 
core. In other words, the non-condensate part can play an important role in 
the core region, and the real particle density distribution can deviate from 
the solution of the GP equation significantly]^ 



2. PROJECTION METHOD 

In order to properly address the issue on quantum fluctuations, I use 
the projection method to construct a many-body wave function for a vor- 
tex by spreading out the vorticity over a finite region. This method was 
first proposed by Hill and WheelerE-a in the context of rotational states of 
a nucleus, and was developed in detail by Peierls, Yoccoz, and ThoulessEl 
a few years later. Recently, this method has been successfully applied to 
a study onthe collective motion of avortex in a uniform two-dimensional 
superfluidity In Madison's experiments on a rotating single-component con- 
densate, the aspect ratio of the condensate between the longitudinal size (z 
axis) and the transverse size (x-y plane) is about 20 : 1. Therefore, it is 
reasonable to consider an effective two-dimensional theory, which assumes 
the structure of a static vortex is invariant under translations along the lon- 
gitudinal direction. The Kelvin wavesE-0 propagating along a vortex filament 
are completely outside the scope of the present study. 

In the mean-field theory, a singly quantized vortex line located at a 
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point ro near the z axis is described by a Feynman wave function]} 



tfp(ri,...rA;r 
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i=i 



^ gs ( ri ,...r^) , (1) 



where N is the total number of atoms in this effective two-dimensional sys- 
tem, \I/g S is the ground state wave function, g(\r — ro|) is solved by the GP 
equation, and 0(r; ro) is a scalar potential describing the circulating velocity 
field. For all practical purposes, the normalized ground state wave function 
is given by the Thomas-Fermi approximation, 
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where R is the radius of the condensate. The normal component of the 
velocity field has to vanish at the boundary, and the velocity potential can 
be easily solved as 

by using an image vortex with an opposite circulation located at (i?/ro) 2 ro. 

As I already mentioned in the introduction, the difficulty of the mean- 
field theory was that the radial wave function g(r) had to vanish at the 
vortex axis because the phase of the single-particle wave function was not 
well-defined. To go beyond the mean-field approximation, I set g(r) to unity, 
and construct an alternative vortex wave function as a linear combination 
of Feynman wave functions centered at different positions, 

*(ri,...r#) = |d 2 ro/(ro)e^ ^ F (ri,...r^;r ), (4) 

where the weighting function /(ro) can be determined from a variational 
calculation.0 This wave function \& has several salient features. First of 
all, it is an eigenfunction of the total angular momentum operator L z with 
eigenvalue Nh. Secondly, it is consistent with the symmetry of the Hamilto- 
nian that the wave function ^f, unlike the Feynman wave function typ, does 
not explicitly contain the vortex position ro. Thirdly, quantum fluctuations 
have been properly addressed, and the single-particle aspect of the vortex 
motion can be described through the weighting function /(ro). 

To compute the weighting function, consider the variational integral, 

= = /d 2 r /d 2 r /(ro)/(r )e^o^ F K)|g|^F(r )) 

/d 2 r /d 2 ro/(r )/(ro) e ™^ F (r )|M/ F (r )) ' 
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where A6q = #o — $o and H is the Hamiltonian. The requirement that (E) 
be a minimum leads to an integral equation for f(ro), 

dr Kp ,r )f(r ) - (E) [ dr J(r' ,r )f(r ) = , (6) 

JO 

where 



o Jo 



p2tt r2/T 

J(r' ,r ) = r dO'J d6 e iN ^ {^Wo)\^{ro)) , (7) 
Jo Jo 

K(r' ,r ) = r d9'J d9 e iNAe ° (n(r' )\H\^ F (r )} . (8) 
Jo Jo 

Further mathematical details on solving this integral equation is given in 



Ref. 16. I will only show the relevant results here. The most important 
quantity in calculating the structure is the overlap between two Feynman 
wave functions. In the limit of large number of particles, the overlap takes 
the following form, 



e™ (^(r[ ) )|^ F (r ))=exp 



1 2R \d 2 r xrj 
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, (9) 



where d = \tq — r' Q \ is the separation between the two Feynman wave func- 
tions, a = 0.424 is a numerical constant, and a = (R 2 /2N) 1 / 2 is approxi- 
mately the two-dimensional interatomic spacing of the ground state in the 
center of the gas. The imaginary part in the overlap is related to the Mag- 
nus force, and makes the weighting function localized in space. For the case 
of a stationary vortex along the z axis, the weighting function is given by a 
simple Gaussian /(ro) ~ exp(— 7q/2<7 2 ) apart from a normalization constant. 



DENSITY AND VELOCITY PROFILES 



Using the wave function in Eq. (Q), one can calculate the number density 
and the current density. The number density is defined as 
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Since the density is independent of the azimuthal angle by symmetry, I can 
perform an angular average, and rewrite the integrand as a function of A#o- 
The radial profile of the number density is 
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x[o] ~ r[o] 



Fig. 1. The left panel shows the column density profile of a vortex in the 
center of the trap. The central density drops about 47% from its ground 
state value. The right panel shows the velocity distribution. 



where 



5(r ,r' ) = /(ro)/(r^ 



JNA6 



(*Hr'o)|*F(ro)} . 



(12) 



The angular average integral in En. (11) can be worked out explicitly as 
a combination of elliptic integrals.!!^ The rest integrals can then be carried 
out numerically. In order to compare with the experimental result I use 
the parameters in Ref. [?]. The scattering length for 87 Rb is 5.45 nmJ!3 The 
frequency of the harmonic trap in the transverse direction is wj_/27r = 219 
Hz, and in the longitudinal direction is u z /2tt = 11.7 Hz. For 10 5 atoms, 
the maximum transverse size is R = 2.6 /xm, and the central density for 
the ground state is about n(0) = 177/im -3 . However, one has to find the 
effective length scale a in the two-dimensional theory from the parameters in 
three dimensions. Since a determines the size of the vortex core, a reasonable 
estimation is the healing length £. For a uniform gas, the healing length is 
given by £ = (87rna) -1 / 2 , which is about 0.2/xm if the central density n(0) 
is used. In units of this effective length scale, a ~ £, the transverse size 
of the condensate is about R/a ~ 13. I have carried out a calculation 
for R/a = 15, and the normalized density profile is shown in Fig. [p. The 
radial dependence of the condensate wave function has been changed from 
the pure two-dimensional profile (1 — r 2 /R 2 ) to the column density of a 
three-dimensional condensate (1 — r 2 / R 2 )^l 2 . 

One can compute the current density in a similar way. The azimuthal 
current density is defined as 



Mr) 



N 
3=1 
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(13) 



The radial distribution of the azimuthal velocity, vg(r) = Jo{r) / p(r), is also 
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shown in Fig. || together with the velocity profile h/mr of a singular vor- 
tex line. In conclusion, a many-body wave function incorporating quantum 
fluctuations is presented to describe a vortex in a BEC. The calculated par- 
ticle density at the vortex axis is finite, and agrees quantitatively with the 
experimental observation. The finite density at the vortex axis cannot be 
described by the GP equation, and manifests the quantum depletion. 
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